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I.  INTRODUCTION 


Since  1957  the  continuous  use  of  space  has  resulted  in  the  buildup  of  a 
large  number  of  space  objects  such  as  active  and  spent  payloads,  rocket 
bodies,  and  miscellaneous  debris,  including  numerous  explosion  fragments.  A 
catalog  of  an  estimated  radar- trackable  population  of  some  5000  satellites  is 
maintained  by  the  North  American  Air  Defense  Command  (NORAD)  and  other 
agencies.  These  objects  are  a  subset  of  an  unknown  but  larger  population  of 
objects  including  those  that  are  too  small  to  be  tracked  by  radar,  but  large 
enough  to  cause  damage  to  another  spacecraft  by  collision. 

In  view  of  this,  much  interest  and  concern  have  been  expressed  recently 
by  several  government  agencies  and  international  organizations  regarding  the 
orbiting  debris  problem  and  the  associated  collision  hazard.  Questions  on  the 
subject  have  been  posed  by  the  United  Na.'ions,  NASA  and  the  USAF/SD  which  have 
resulted  in  some  assessments  of  the  collision  hazard  at  present  and  in  the 
near  future.  Also,  work  is  currently  underway  to  revise  the  USAF/SD 
Regulation  550-11  to  establish  Space  Division  guidelines  and  responsibilities 
for  program  directors  and  offices  involved  in  the  management  of  geosynchronous 
satellites. 

The  purpose  of  the  present  study  is  to  update  the  existing  data  on  the 
orbiting  population  (Ref.  1),  and  to  reexamine  the  collision  hazard  issues  as 
currently  perceived.  Probabilities  of  collision  at  low  altitudes  and  in  the 
geosynchronous  corridor  will  be  estimated  using  approximate  methods  developed 
in  the  study.  Potential  operational  and  design  approaches  that  can  aid  in 
reducing  the  space  debris  population  growth  rates  will  be  outlined.  It  is 
hoped  that  these  efforts  will  contribute  towards  a  better  understanding  and 
awareness  of  the  collision  hazard  so  that  appropriate  action  can  be  taken  to 
maintain  a  relatively  low  risk  environment  for  current  and  future  satellite 
systems. 
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II.  SPACE  OBJECT  POPULATION 


A.  GROWTH  RATE 

On  the  basis  of  a  worldwide  network  of  sensors,  NORAD  tracks  and 
maintains  a  catalog  of  most  of  the  manmade  orbiting  objects  in  space. 

Although  NORAD  is  the  primary  source  of  orbit  elements  for  these  objects, 
other  sources  are  NASA  Satellite  Situation  Reports,  RAE  tables  of  earth 
satellites,  TRW  Space  Log  and  various  military  and  civilian  agencies  including 
those  of  the  United  Nations. 

The  orbiting  population  of  objects  consists  of  operational  and  decayed 

payloads,  miscellaneous  mission-related  debris  such  as  rocket  bodies,  clamps, 

shrouds,  etc. ,  and  a  multitude  of  explosion  fragments.  The  latter  constitute 

more  than  60%  of  the  trackable  population  (Ref.  2).  For  example,  the  number 

and  type  of  objects  listed  in  the  NORAD  catalog  for  27  April  1980  are  given  in 

Table  1  (Ref.  3).  It  should  be  noted,  however,  that  the  cataloged  population 

includes  only  those  objects  which  can  be  tracked  by  radar  (greater  than  0.01 
2 

m  in  area)  up  to  5000  km  in  altitude  or  optically.  The  currently  estimated 

2 

total  population  (down  to  1  cm  in  size)  is  believed  to  be  on  the  order  of 

10,000  to  15,000  objects.  In  the  geosynchronous  corridor,  the  total 

2 

population  down  to  1  cm  is  probably  as  much  as  an  order  of  magnitude  greater 
than  that  which  is  cataloged.  More  detail  on  the  synchronous  satellite 
catalog  is  provided  in  the  Appendix. 

The  approximate  growth  rate  of  the  cataloged  population  is  exhibited  in 
Fig.  1.  The  observed  rate  of  increase  (on  the  basis  of  the  smoothed  data  of 
Fig.  1)  is  between  9  and  13%  per  year  depending  on  the  time  period  of 
interest.  An  instructive  computer  illustration  of  the  space  objects  which 
could  be  seen  by  an  observer  in  a  300  nmi  circular  orbit  for  a  period  of  10 
seconds  is  displayed  in  Fig.  2  where  the  apparent  size  of  the  object  (circle 
radius)  denotes  relative  distance  from  the  observer. 
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Table  1.  NORAD  Catalog  of  All  Objects  for  27  April  80 


Earth  Orbit  Payloads 

1055 

Earth  Orbit  Debris 

3384 

Deep  Space  Payloads 

61 

Deep  Space  Debris 

52 

Current 

Objects 

4552 

Decayed  Payloads 

1399 

Decayed  Debris 

5831 

Decayed 

Objects 

7230 

Total  Objects 

Synchronous  Population  -  28  April  80 


Satellites 

DOD(SCF)/NATO  23 
NASA  1 7 
COMSAT  CORP  1  3 
LINCOLN  LABS,  RCA,  W.U.,  ESA  10 
USSR  16 
UK,  CAN,  JAP,  FR,  IT,  INDONESIA  18 


97 

Rocket  Bodies  and  Old  Satellites 


With  Current  Tracking 

47 

No  Recent  Tracking 

56 

103 

Total  Synchronous 


« 


4552 

7230 

11782 


97 


j_03 

200 
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REVISED  TOTAL  OF  ALL  OBJECTS 
LARGER  THAN  8  cm  DIA.  (SOURCE: 
NASA  JSC  BRIEFING  MARCH  1979! 


1970 


TOTAL  TRACKED  OBJECTS 
(NORAD  CATALOG) 


/ 

/ 

/ 


/ 


/ 


GEOSYNCHRONOUS  ORBIT 
TRACKED  OBJECTS 


1975  1980  1985  1990 


Cataloged  Objects  Population  History 


Fig.  2.  Computer  View  of  a  Sample  of  Space  Objects 
in  a  300  nmi  Circular  Orbit 


B. 


ORIGIN  AND  DISTRIBUTION 


The  principal  sources  of  the  orbiting  population  of  objects  to  date  have 
been  the  nearly  2000  launches  by  the  U.S.,  the  USSR,  and  several  other 
countries.  The  history  of  known  launches  and  payloads  placed  in  orbit  since 
1975  is  exhibited  in  Fig.  3.  The  principal  U.S.  launch  vehicles  are  the 
Delta,  Atlas/Agena,  Centaur,  Titan  III  and  the  Scout  stages  which  have 
accounted  for  the  estimated  33  U.S.  launches  in  1978.  The  total  of  87  USSR 
launches  in  1978  have  resulted  in  a  large  number  of  the  COSMOS  series  of 
satellites  some  of  which  were  returned  to  earth  or  commanded  into  a 
destructive  reentry  into  the  atmosphere.  By  1  January  1979  there  were  six 
launches  bv  France,  four  by  Japan,  two  by  China  and  one  by  the  U.K.  From  1988 
the  vearlv  total  number  of  launches  varied  between  108  and  128  for  all 
countries  (Ref.  A). 

Payload  and  rocket  body  explosions  in  orbit  are  another  source  of  space 
dehris.  It  is  known  for  example,  that  there  have  been  over  50  such  explosions 
as  of  4  October  1978.  A  number  of  spacecraft  explosions  are  known  to  have 
taken  place  in  orbit  related  to  antisatellite  tests  (ASATs)  and/or  attempts  to 
inject  payloads  into  geosynchronous  orbit.  A  recent  example  of  the  latter  was 
the  apparent  loss  of  the  RCA  Satcom  3  and  the  Japanese  ECS  2  satellites  which 
seem  to  have  experienced  anomalous  burns  of  the  apogee  injection  motors  that 
resulted  in  a  loss  of  telemetry  before  the  completion  of  the  insertion  burn. 
Another  example  of  such  a  malfunction  is  the  SKYNET  IB  launch  on  19  August 
1970  on  a  Thor-Delta  vehicle  which  placed  a  communication  payload  into  a  137  x 
20099  nmt  transfer  orbit.  On  22  August  1970,  the  apogee  burn  motor  (ABM)  was 
fired  which  should  have  placed  the  satellite  into  a  nearlv  synchronous 
equatorial  orbit.  Approximately  14  seconds  after  the  initiation  of  the  apogee 
hum,  all  tracking  and  telemetry  data  were  abruptly  lost.  Subsequent  analyses 
of  the  data  (Refs.  5  and  8)  indicated  that  the  ABM  performed  anomalously 
resulting  in  "regressive  burning."  An  estimated  AV  =  4000  fps  out  of  the 
required  AV  =  5814  fps  appears  to  have  been  added  at  apogee.  The  resultant 
orbit  for  the  assumed  explosion  fragments  was  estimated  as  approximately 
5192  x  20088  nmi  inclined  5.30°  to  the  equator.  The  ve 1 oc i t v  at  apogee 
relative  to  the  geosynchronous  velocity  is  about  1700  fps  which  mav  represent 


a, 


NUMBER,  LAUNCHES  AND  PAYLOADS 


a  significant  collision  hazard  to  any  geosynchronous  object  in  the  vicinity  of 
this  orbit.  The  number  and  distribution  of  the  explosion  fragments,  if  an 
explosion  in  fact  occurred,  cannot  be  easily  determined.  However,  the  high 
perigee  altitude  of  the  final  orbit  precludes  any  significant  reduction  of  the 
debris  because  of  atmospheric  or  other  perturbative  effects.  Figure  4 
exhibits  the  orbit  geometry  for  this  event. 

Other  potential  sources  of  space  debris  are  collisions  between  debris 
objects  and  spacecraft  or  other  debris  objects.  It  is  considered  likely,  for 
example,  that  the  solar  panel  malfunction  on  the  geosynchronous  GEOS-2 
spacecraft  is  a  result  of  a  collision  with  a  small  object  (Ref.  7).  Also,  the 
breakup  of  PAGEOS  balloon  in  July,  1975  after  9  years  in  orbit  appears  to  have 
resulted  from  a  similar  cause.  The  low  orbit  COSMOS  954  was  similarly 
reported  to  have  suffered  a  sudden  depressurization  thought  to  be  caused  by  a 
collision  with  debris  (Ref.  2).  Air  Force  satellite  programs  in  synchronous 
orbit  have  reported,  on  two  separate  occasions,  sudden  small  changes  in 
spacecraft  angular  momentum.  These  unexplained  changes  may  have  been  caused 
by  impacts  with  natural  or  manmade  objects. 

The  possibility  of  a  self-perpetuating  debris  belt  caused  by 
intercollisions  between  objects  in  space  also  exists.  This  process  would 
parallel  certain  theories  concerning  the  growth  of  the  asteroid  belt.  The 
debris  flux  in  such  a  belt  could  exceed  the  natural  meteoroid  flux  in  the  not 
too  distant  future  if  present  trends  continue  (Ref.  8). 

C .  DENSITY  OF  OBJECTS  AS  A  FUNCTION  OF  ALTITUDE 

1.  LOW  ALTITUDE  ORBITS 

A  total  of  4174  space  objects  from  the  NORAD  Catalog  for  April  28,  1980 
have  been  examined  numerically  to  determine  the  distribution  with  altitude. 

The  results  are  presented  in  Figs.  5  and  6.  It  can  be  seen  that  the  majority 
of  satellites  are  within  1.5  earth  radii  (E.R.)  with  a  small  increase  to 
synchronous  altitude  where  an  increment  occurs  at  (5.62  E  R.  Figure  6 
indicates  that  the  greatest  density  of  objects  occurs  between  500  and  1500  km 
altitude.  Distributions  by  semimajor  axis  A  (normalized  with  respect  to  earth 
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Fig.  4.  SKYNET  IB  Geosynchronous  Orbit  Injection  Geometry 
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Cumulative  Number  of  Objects  va.  Geocentric  Radius  In 
Units  of  Earth  Radius 


Density  vs.  Altitude 


radius),  eccentricity  E,  ascending  node  (NODE)  and  the  argument  of  perigee 
(OMEGA)  are  given  in  Fig.  7.  It  can  he  seen,  for  example,  that  there  are  3791 
obiects  with  A  greater  than  1  hut  less  than  2,  and  F.  greater  than  zero  hut 
less  than  0.01.  The  distribution  with  respect  to  orbital  inclination  (INC)  is 
given  in  Fig.  8. 

2.  GEOSYNCHRONOUS  ORBITS 

The  geosynchronous  region  of  space  is  a  unique  corridor  where  most  of  the 
present  and  future  communication  satellites  are  or  will  he  located.  The 
principal  reason  for  this  is  that  these  satellites  appear  to  be  either 
stationary  with  respect  to  the  earth  or  are  slowly  drifting  relative  to  some 
initial  longitude.  The  satellites  are  stationary  when  thev  are  in  circular 
orbits  in  the  equatorial  plane  with  periods  of  revolution  equal  to  the  rota¬ 
tional  period  of  the  earth.  Secular  or  continuous  drifting  occurs  due  to 
slight  differences  in  the  orbital  periods  of  revolution.  Periodic  drifting 
may  also  occur  because  of  orbital  eccentricity  or  inclination  effects  as  well 
as  perturbing  effects  of  the  geopotential  field.  For  example,  the  J2  2 
tesseral  harmonic  causes  an  oscillatory  (1 i brat i onal )  motion  about  the  geo- 
potential  stable  points  at  75.3°E  and  285. 3°F  longitudes  with  long  periods  (on 
the  order  of  years)  of  lihration  (see  Fig.  9).  A  large  number  of  obiects 
lihrating  about  these  points  mav  pose  an  additional  hazard  to  the  active 
payloads  placed  at  or  near  these  locations.  The  longitudinal  location  for  134 
geosynchronous  obiects  listed  in  the  appendix  is  illustrated  in  Fig.  10. 

NORAD  publishes  the  orbital  elements  on  these  obiects  on  the  basis  of  data 
obtained  from  NORAD,  the  Satellite  Control  Facilitv  (SCF),  NASA  and  Tntelsat 
tracking  networks.  The  data  on  the  orbital  elements  is  provided  hv  NORAD  in 
the  form  of  punched  cards  received  via  TWX  on  a  dailv  basis.  Updates  for  anv 
given  object  are  provided  typically  within  a  one  or  two  week  period. 

Figure  ll  presents  the  longitudinal  location  of  the  current  and  planned 
communication  satellites  in  geosynchronous  orbit  (Refs.  9  to  11).  Manv  of 
these  comsats  are  station-kept  within  small  longitudinal  and  latitude  hands 
while  others  drift  as  the  result  of  natural  forces.  A  few  havp  also  been 
removed  from  geosynchronous  orbit  after  the  completion  of  their  mission  (e.g. , 
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raction  of 


Geopotent 


Intelsat  3,  F-2,  F-4  and  F-6  to  altitudes  of  400  to  3700  km  above  synchronous 
a  It  It  ode )(Ref .  12).  The  ATS-F  has  been  removed  from  the  synchronous  altitude 
orbit  to  an  altitude  250  nmi  below  synchronous.  The  distribution  of  a  sample 
of  134  geosynchronous  objects  with  inclination  and  eccentricity  is  given  in 
Figs.  12  and  13,  respectively.  These  results  indicate  that  typical  geo¬ 
synchronous  orbits  are  nearly  circular  with  low  inclinations  relative  to  the 
equator.  A  few  satellites  (e.g.  ELECTRON,  and  their  rocket  bodies  have 
inclinations  above  60°  (see  Appendix).  Object  densitv,  illustrated  in  Fig. 

14,  decreases  as  the  altitude  from  the  geosynchronous  orbit  and  latitude 
increase.  The  results  in  Fig.  14  were  obtained  numerically  by  taking 
"snapshots"  of  objects  once  each  hour  over  a  period  of  twenty-four  hours 
within  a  torus  of  latitude  band  A<f>  and  altitude  band  Ah  centered  about  the 
geosynchronous  altitude.  The  total  number  of  observations  was  then  summed  and 
divided  by  24  to  get  an  average  number  of  sightings  in  each  band.  The  average 
density  is  the  number  of  sightings  divided  by  the  volume  of  the  torus.  The 
nonuniformity  of  the  inclinations  of  the  objects  orbits  within  each  latitude 
band  is  illustrated  in  Fig.  14  by  the  average  inclination  i  as  noted. 

Tiie  ascending  node  distribution  for  154  geosynchronous  objects  presented 
in  Fig.  15  indicates  that  the  orbital  planes  are  not  distributed  uniformly, 
but  tend  to  be  concentrated  in  certain  regions  of  space.  The  impact  of  this 
is  to  increase  the  probability  of  collision  for  some  satellites  above  t he* 
representative  figures  calculated  for  the  synchronous  class  of  satellites. 
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Number  vs.  F.ccentri  ci  tv  for  Geosynchronous  Obiects 


A<|>  =  LATITUDE  BAND,  deg 
I  -  AVERAGE  INCLINATION  IN  A <f> 


Ah,  DISTANCE  FROM  GEOSYNCHRONOUS  ALTITUDE,  nmi 


Ftg.  14.  Ob.iect  Density  in  Geosynchronous  Orbits 
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Right  Ascension  of  the  Ascending  Node  in  Geosynchronous  Orbits 


III.  ANALYSIS  AND  RESULTS 

A.  UNIFORM  DENSITY  METHOD 

The  probability  that  any  two  satellites  will  collide  is  generally  a 
function  of  their  orbital  parameters,  size  and  time.  The  collision  cannot 
take  place  until  the  orbits  intersect  or  approach  each  other  to  within  the 
dimensions  of  the  satellites  of  interest.  This  may  occur  even  for  initially 
nonintersecting  orbits  because  of  the  effects  of  the  earth's  oblateness,  air 
drag  and  solar-lunar  perturbations  which  alter  the  orbital  parameters  in  time. 

An  approximate,  yet  reasonably  accurate  method  for  estimating  the 
probability  of  collision  between  a  target  satellite  of  projected  area  Ac  and  a 
debris  object  is  based  on  the  assumption  that  the  space  object  density  p  is 
uniform  (constant)  and  that  the  target  area  Ac  sweeps  out  a  volume  AcAt  in 
time  At.  It  can  be  shown  (Ref.  1)  that  the  fractional  number  of  objects 
encountered  in  the  volume  A£  At  is  then  equal  to  the  probability  of  collision 
if  it  is  much  less  than  unity.  Thus, 

p(col)  =  pv  A  At  ( 1 ) 

r  c 


where  vr  is  the  average  relative  velocity  between  the  satellite  and  the  set  of 
objects  considered.  For  example,  vr  7  km/sec  for  circular  orbits  up  to 
2000  km  in  altitude  (Ref.  8). 


The  relative  velocity  vr  at  the  trace  intersection  between  any  two 
satellites  in  orbit  planes  with  a  mutual  inclination  of  Ai  can  he  expressed  as 
(Ref.  1): 
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Here 


where  rx  is  the  radius  at  trace  intersection  and  Aj  =  aj/rx  ^  -  a2^ry.' 
ap  >  ei>  e2>  h  are  the  semimaior  axes,  eccentricities  and  the  gravitational 
constant,  respectively. 

For  the  case  of  a  circular  and  an  eccentric  orbit  (e.g.  ~  0,  rx  =  82) 

Eq.  (2)  becomes 


(3) 


For  circular  orbits  of  equal  period  (i.e.,  aj  =  a2>  ei  =  e2  =  °) »  Kq .  (2) 
reduces  to 

.  Ai 

v  =  2v  sin  — 
rc  c  2 

(A) 

“v  sinAi  for  Ai  small, 

c 

which  is  equal  approximately  to  the  cross-track  (normal)  component  of  relative 
velocity  at  encounter  with  vc  =  circular  orbital  velocity.  An  average 
relative  velocity  vrc  for  this  case  can  be  defined  as  an  average  value  of  vrc 
over  a  mean  anomaly  range  — tt / 2  <  M  <  it/2,  i.e. 


v  K  —  sin  Ai  I  cos  MdM 
rC  11  J-  tt/2 

2v 

"  — -  sin  Ai 

TT 


(5) 


U 


An  approach  which  can  be  applied  to  the  case  of  a  spacecraft  in  a 
circular,  inclined  orbit  is  based  on  the  determination  of  obiect  density  in  a 
spheroidal  torus  containing  the  spacecraft  orbit  (Fig.  16).  This  approach  is 
similar  to  that  used  in  the  theory  of  interplanetary  encounters  (Ref.  13). 

The  spheroidal  torus  is  defined  by  the  relative  orbit  plane  inclination  Ai  and 
the  spacecraft  orbit  radius  R.  The  probability  of  the  object  colliding  with 
yhe  spacecraft  in  a  time  interval  At  is 


p(col/At) 


.  -AAt 
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«  AAt 


for  AAt  <<  1 
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where 
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p  v  A 
r  c 
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T 


r _ 

sin  Ai 


=  impact  rate  per  unit  time 
T  =  obiect  orbit  period  of  revolution 

S  =  Rs/R 

Rg  =  spacecraft  radius 

vr  =  relative  velocity  at  encounter 

I ur I  =  absolute  value  of  the  radial  component  of  vr 
For  a  circular  target  orbit  and  an  eccentric  object  orbit  (Refs.  1,  13) 
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(VOLUME  =  4ttR2  dR  sin  Ai) 


Fig.  11).  Toroidal  Geometry  for  Geosynchronous  Orbits 


where  Also,  ap  a-^  >  are  c^e  semimajor  axes  of  the  object  and 

target  orbits,  and  e^  is  the  eccentricity  of  the  object  orbit. 

The  probability  of  the  target  spacecraft  colliding  with  N  objects  in  the 
torus  (defined  by  an  average  inclination  Aigv  of  the  objects'  orbits  relative 
to  that  of  the  target  orbit)  is  equal  to  N  times  that  given  by  Eq .  (6).  For 
example,  if  a  set  of  62  geosynchronous  objects  is  considered  with  an  average 
eccentricity  ej  =  0.00113  and  orbit  plane  inclination  Ai  =  1.07°  relative  to 
the  equatorial  plane,  then  the  velocity  relative  to  a  spacecraft  in  a  circular 
geosynchronous  equatorial  orbit  at  encounter  is  from  Eq.  (3)  vr  =  57.5  m/sec 
where  a^  =  1  and  rx  =  geosynchronous  radius.  Consequently,  a  1000  day 
probability  of  collision  for  a  spacecraft  of  50  m  radius  is  by  Eq.  (6)  equal 
to  3.88  x  10“5. 

Using  Eq.  (1),  vr  =  7  km/sec,  which  is  the  average  relative  encounter 
velocity  between  satellites  below  2000  km  (Ref.  8)  and  the  spatial  density  p 
in  Fig.  7,  the  probability  of  collision  for  circular  low  altitude  missions  of 
1000  days  is  given  in  Fig.  17  for  the  trackable  population  of  4174  objects  and 
an  estimated  population  of  8400  objects  of  smaller  cross-section.  The  results 
show  that  the  1000  day  probability  of  collision  for  a  very  large  spacecraft 
(~50  m  radius)  is  on  the  order  of  4  to  8%  and  for  a  20  m  radius  satellite  it 
is  on  the  order  of  0.6  to  1.2%.  However, for  smaller  spacecraft  (~3  m  radius), 
current  probability  of  collision  is  an  order  of  magnitude  lower.  The  results 
for  a  sample  of  62  to  620  objects  in  the  geosynchr onous  corridor  are  given  in 
Fig.  18  where  an  average  relative  velocity  of  only  36.5  m/sec  was  assumed. 

This  velocity  is  the  north-south  (inclination)  component  of  the  average 
relative  velocity.  Average  collision  probabilities  in  the  geosynchronous 
orbit  are  thus  seen  to  be  about  three  orders  of  magnitude  lower  than  those  at 
low  altitudes.  However,  the  effects  of  density  variation  with  altitude, 
inclination  (latitude)  or  longitude  are  not  included  in  these  results. 

B.  VARIABLE  DENSITY  METHOD 

Figure  14  shows  that  the  observed  object  density  p  ( h ,  )  is  a  function  of 
height  h  above  or  below  geosynchronous  orbit  and  latitude  <f>.  Since  a 
satellite  whose  orbit  plane  is  inclined  to  the  equatorial  plane  spends  only  a 
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Fig.  17.  Probability  of  Collision  -  Low  Earth  Orbits 
(800  to  1500  km)  for  1000  Days 
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Fig.  18.  Probability  of  Collist  on  -  Geosynchronous 
Orbits  for  1000  Days 
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fraction  of  its  period  of  revolution  T  near  the  equator  where  the  density  is 
greatest,  the  probability  of  collision  is  similarly  affected.  In  general,  the 
probability  of  collision  in  any  altitude  band  Ah  =  1^  -  hj,  per  one  revolution 
of  the  satellite,  is 


p(col/rev) 


A  v  T 
c  r 


p(h ,  4>)f  ( )d<J>dh 


(8) 


where  p  (h ,  )  is  the  object  density  function,  f((J>)  is  a  weighting  function 
which  can  be  derived  from  the  time  fraction  spent  by  a  satellite  in  a  latitude 
band  Aj>  as  shown  in  Fig.  19  and  vr  is  the  average  velocity  of  the  target 
satellite  relative  to  the  objects  of  interest.  For  the  case  of  a  target 
satellite  with  inclination  i  and  a  set  of  geosynchronous  objects  with  an 
average  inclination  i  and  average  right  ascension  of  the  ascending  node  ft  , 
the  angle  Ai  between  the  target  orbit  plane  and  that  of  the  "average  plane"  of 
the  objects  is 


Ai  =  cos  (cosi  cosi  +  sini  sini  cosAft) 


(9) 


where 

Aft  =  I  ft  -  ft  I 

ft  =  right  ascension  of  the  ascending  node  of  the  target  satellite. 

Since  Ai  will  in  general  vary  in  time  because  of  external  perturbations, 
the  average  value  of  Ai  is 


Ai 

av 
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(10) 
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DWELL  FRACTION  t/T,  PERCEN 


i  -  SATELLITE  ORBIT 
INCLINATION 


T  =  ORBIT  PERIOD 
t  =  TIME  IN  LATITUDE 
BAND  ±<f> 


Hig.  19.  Dwell  Fraction  vs.  Latitude  For 
Circular  Orbits 


where 


Ai  =  i  +  i  when  Aft  =  it  (11) 

max 


and 


At  .  =  1 1  —  i  I  when  Aft  =  0  (12) 

mm 


Therefore , 


Ai 


av 


j  i  for  i  >  i 
I  i  for  i  <  i 


(13) 


The  average  relative  velocity  vr  for  the  case  of  nearly  circular 
geosynchronous  orbits  can  be  approximated  by  Eq .  (5)  because  the  normal  (nut- 
of  plane)  or  north-south  component  predominates. 

Equation  (8)  has  been  evaluated  approximately  for  circular  orbits  with 
i  =  1°,  3°  and  >°  and  object  densities  of  Fig.  14  for  different  altitude  bands 
above  and  below  geosynchronous  orbit.  The  results,  given  in  Fig.  20,  are 
valid  for  a  typical  small  satellite  (effective  collision  radius  =  20  ft 
(b.l  m))  and  a  mission  duration  of  1000  days.  The  results  for  other  values  of 
Rcj  are  proportional  to  the  square  of  the  effective  radius  of  collision  and 
time.  The  collision  probability  for  a  target  satellite  in  an  elliptic 
geosynchronous  orbit  can  be  obtained  by  averaging  the  results  of  Fig.  18  as 
shown  in  Fig.  21. 
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Fi£.  20.  Collision  Probability  for  Ohieets  in  Circular 
Geosynchronous  Orbits  in  1O00  Davs 
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DISTANCE  OF  CLOSEST  APPROACH  METHOD 


1.  GENERAL  CONSIDERATIONS 

A  collision  between  two  objects  can  take  place  at  or  near  a  trace 
intersection  point  if  the  distance  of  closest  approach  R^jn  Is  e9ual  to  or  is 
less  than  the  effective  collision  radius  Rg.  For  satellites  in  circular 
orbits  with  a  mutual  orbit  plane  inclination  a,  R^jj^j  occurs  twice  per 
revolution  in  the  vicinity  of  the  intersection  between  the  orbit  planes  (nodal 
axis)  as  illustrated  in  Fig.  22.  For  each  instant  when  satellite  1  is  at  or 
near  the  nodal  axis,  satellite  2  is  at  a  position  4,  3,  2  or  1  corresponding 
to  an  R^f^  for  that  pass.  The  angular  increment  Au  =  nj  I  (T2  -  Ti)l  where  n-j 
is  the  mean  motion  of  satellite  1  and  T^,  T7  are  the  periods  of  revolution  for 
satellites  1  and  2,  respectively.  The  angular  change  Au  per  revolution  is 
typically  a  fraction  of  a  degree  for  geosynchronous  satellites.  For  example, 
^MIN  f7°r  norttl  and  south  bound  passes  of  OPS  6391  (SDC  object  No.  10669)  and 
WESTAR-A  (SCD  object  No.  7250),  is  given  in  Table  2  for  several  days  in  April 
of  1980  as  was  determined  approximately  via  the  numerical  simulation  described 
in  Ref.  14.  Different  R^^N  values  can  thus  be  seen  to  have  occurred  before 
and  after  the  lowest  value  RMIN  =  5.28  nmi  of  21  April  1980.  The  probability 
of  collision  between  the  objects  for  each  pass  is,  in  general,  a  function  of 
RmI N »  t*ie  tracking  uncertainty  a  in  the  location  of  each  object  and  the 
effective  collision  radius  Rg.  The  latter  can  be  defined  as  one  half  the  sum 
of  the  maximum  dimensions  for  both  objects.  For  0  =  0  and  Rvijn  ^  no 

collision  can  take  place.  For  a  f=  0,  there  is  a  nonzero  probabi  1  i  tv  of 
col lision. 

2.  MISS  DISTANCE  AND  POSITION  UNCERTAINTY 

An  approximate  collision  probability  method,  described  brieflv  in  (Ref. 
15),  considers  the  effects  of  the  uncertainties  associated  with  the  three 
dimensions  (coordinates)  of  the  miss  distance  Ry^N’  assumed  that  the 

uncertainties  are  Gaussian  (normal)  with  zero  biases  and  equal  variance  and 
that  they  are  uncorrelated.  The  assumption  is  applied  to  the  position  data  of 
each  of  the  tracked  satellites.  In  view  of  this  and  the  fact  that  the 
orientation  of  the  coordinate  system  containing  the  distance  of  closest 


Table  2.  Close  Encounters  Between  OPS  6391  and  WESTAR-A  Geosynchronous  Satellites 
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approach  at  encounter  is  arbitrary,  a  plane  xy  which  is  normal  to  the  relative 
velocity  vector  and  which  contains  the  vector  of  closest  approach  R^j-^  can 
always  be  found.  The  bi-variate  normal  density  function  for  this  case  can  be 
expressed  as  (Ref.  16) 


f(x,y) 


where 


RMIN 


(X  .  +  Y  .  ) 

'  min  min' 


1/2 


Since  a  collision  can  occur  only  if 


(14) 


min 

X  . 
min 


R  «  y  <  R  +  Y  ,  ) 

s  s  min  I 

R  <  x  <  R  +  X  .  ( 

s  s  min  1 


(15) 


where  X^^,  Ymjn  are  the  coordinates  of  R^j^  and  where  Rs  is  the  effective 
collision  radius  for  both  satellites,  the  probability  of  collision  is 


p(col)  = 


(16) 


f (x,y )dydx 


(17) 


for  Rs  «  Rmin. 

Equation  (17)  represents  the  probability  of  collision  between  two 
satellites  per  encounter  when  the  distance  of  closest  approach  R^j^,  the  1  o 
tracking  uncertainty  in  the  position  of  each  satellite,  and  the  effective 
radius  of  collision  Rg  are  given. 

A  similar  result  can  be  obtained  by  orienting  the  x,  v  coordinates  along 
RMIN  an<^  normal  to  it  respectively.  The  probability  of  collision  is  then 


p(col)  =  pxpv 


(18) 


where 


=  p[  A  <  x  <  B  ] 


=  p[£  <  x 


< 

o 1 


1 

7?tto 


2 

P 


2 


du 


probability  of  collision  along  x 


p [ -Rs  <  y  <  RS1 


p 


V 


Here 


2 

u 


2 


da 


=  probability  of  collision  in  a  direction  normal  to 


A  'VilN  ~  Rs 
B  =  "min  +  Rs 

and  A/a,  B/o  and  R  /o  are  standardized  Gaussian  variables. 

Assuming  now  that  the  probability  of  collision  at  j-th  closest  approach 
is  p.(col)  the  probability  of  miss  is 


p^(miss)  =  1  -  p-(col) 


09) 


The  probability  of  missing  at  all  n  closest  approaches  is 


n 

p(miss/n)  =  n  fl  -  p.(col)l 
1=1  ’ 


■>D 


(20) 


if  all  probabilities  are  independent.  Therefore,  the  probability  that  there 
will  be  at  least  one  collision  during  this  period  of  encounter  is 


p(col)  =  1  -  p(miss/n) 


n 

B  E  p.(col)  (21) 

j-1  J 


when  all  pj(col)  «  1. 

A  plot  of  Eq.  (17)  for  Rg  =  20,  50  and  100  ft  is  given  in  Figs.  22  to  26 
as  a  function  of  o  with  R^^  as  a  parameter.  It  can  be  seen  that  a  maximum 
probability  of  collision  occurs  when 


rmin 
0  =  ~7T 


(22) 


which  is  plotted  as  a  function  of  *n  Fig.  28.  The  use  of  F.qs.  (17) 

and  (21)  is  illustrated  in  Fig.  29  for  the  case  of  the  close  encounters 
between  the  WESTAR-A  and  the  OPS  6391  satellites  given  in  Table  2.  The 
collision  probabilities  for  several  passes  are  plotted  as  a  function  of  o. 

The  maximum  probability  of  collision  for  17  passes  is  given  as  the  sum  of  the 
maximum  probabilities  for  each  pass.  Assuming  no  correlation  between  passes, 
the  maximum  (upper-bound)  probability  of  collision  for  these  encounters  is 
2.0  x  10  6  for  the  8-day  time  period. 


PROBABILITY  PER  PASS 


Collision  ProHnbi 1  i  tv  ns  n  Function  of  Position 
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TV.  SUMMARY  AND  CONCLUSIONS 


A  summary  of  the  present  and  proiected  prohahi 1 i t i es  of  collision  for  a 
ten  meter  radius  spacecraft  in  low  earth  orbit  (LUO)  and  the  geosynchronous 
corridor  (CEO)  is  given  in  Tahte  3.  Tt  can  be  seen  that  the  probability  of 
collision  in  low  earth  orbits  is  currently  on  the  order  of  0.15  to  0.3?  for  a 
1000-dav  mission.  Another  wav  of  interpreting  the  probability  of  0.3?  is  that 
one  of  90  such  spacecraft  would  be  expected  to  experience  a  collision  in  10 
years.  Probabilities  for  the  1 98 S  to  1995  time  frame  are  on  the  order  of  2.5 
to  7.3  times  current  values  on  the  hasis  of  an  assumed  debris  growth  rate  of 
13?  per  vear.  On  the  other  hand,  the  results  for  the  geosynchronous  orbit  are 
some  three  orders  of  magnitude  lower  in  value  mainly  because  of  lower  relative 
velocities  at  encounter.  These  results,  however,  do  not  include  the  effects  of 
ascending  node  concentrations  at  certain  longitudes.  The  calculated  probabil¬ 
ities  of  collision  at  low  altitude  are  also  likely  to  be  too  low  due  to  the 
noninclusion  of  many  smaller  particles  believed  to  he  a  result  of  satellite 
explosions  and  ASAT  tests. 

The  picture  emerging  from  this  study  indicates  that  the  cluttering  of 
space  with  debris  must  be  reduced  in  the  future  to  minimize  the  collision 
hazard  and  improve  spacecraft  survivability.  It  appears  that  the  current 
debris  growth  rates  of  9  to  13%  are  excessive  and  that  the  associated 
collision  hazards,  particularly  for  larger  spacecraft,  will  become 
unacceptable  if  these  growth  rates  remain  unchanged.  There  is  little  doubt 
that  larger  satellites  and  space  stations  will  be  built  in  the  future,  because 
of  increasing  demands  for  more  and  better  spac"  communications,  meteorology, 
navigation,  and  remote  sensing. 

The  question  of  "what  solutions  can  be  postulated"  can  be  answered  onlv 
by  resolution  of  the  space  debris  issues  falling  into  three  categories:  1) 
satellite  and  vehicle  design,  2)  operational  procedures  and  practices,  and  3) 
national  and  international  policies  and  treaties.  Tn  the  area  of  vehicle 
design  the  principal  approaches  should  consider  space  systems  for  litter-f;.-< 


Table  3.  Summary  of  Present  and  Projected  Collision 
Probabilities  for  a  1000-Day  Mission 


in  x>  r>| 

GO  O  6 

§ 

o 


1 

X 

X 

o 

m 

o 

00 

o 

o 

AJ 

o 

4-1 

4-1 

m 

A-> 

X> 

\ 

m 

1 

i 

o 

Csl 

1 

o 

O 

- — i 

o 

■  ' 

X 

X 

m 

X 

• 

4J 

i— > 

•f— 

X 

u 

c 

4-> 

w 

*r— 

3 

O 

V-4 

/"■\ 

✓**s 

C 

o 

X 

Xl 

o 

V-/ 

C/} 

to 

V- 

,n 

o 

3 

O 

o 

3  o 

,  r* 

4-1 

VI 

U3 

■.-i 

Ed 

to  LL5 

•r-»  U3 

O 

U 

o 

TD 

O 

3  -J 

"O  -J 

C 

(0 

tO 

•— < 

to 

>, 

w 

xs 

c 

U 

c 

•u  c 

c 

t/J 

X 

■f4 

•H 

X  -H 

•r— 

o 

5 

u 

U 

»- 

a> 

0 

4-> 

0) 

4J 

AJ 

(U  -u 

o 

^4 

u 

4-1 

U- 

Vj  14_j 

4-J  U— 

QJ 

X 

OJ 

to 

0)  CO 

to 

i 

1 

4-> 

u 

B 

1- 

4-)  U 

G  ^ 

0> 

u 

u 

a;  U 

u 

o 

O 

E 

a> 

V 

G  0) 

>  .  a> 

UJ 

U2 

u 

4J 

o 

u 

4->  u 

o 

-J 

C 

to 

<4-J 

X 

C  to 

u -  to 

a. 

*i“ l 

a . 

<u  a. 

•*-»  a. 

H 

t/> 

Cju 

to 

H  to 

U-  C/) 

/—V 

(0 

-D 

separation,  reusability,  retrlevabl llty,  earth  escape  or  destructive  reentrv 
Into  the  atmosphere.  The  feasibility  of  debris  collection  systems  should  be 
examined. 

In  the  operational  approach  and  procedure  areas,  consideration  should  be 
given  to  satellite  separation  techniques  (use  of  nontntersecting  orbits, 
etc.),  avoidance  of  crowded  regions  of  space,  and  the  disposal  of  spent 
satellites  to  "graveyard"  orbits.  Accidental  or  deliberate  destruction  of 
spacecraft  should  be  minimized.  Some  of  the  possible  approaches  are 
summarized  in  Fig.  30.  In  short.  Planet  Earth  appears  to  need  a  consistent 
and  universally  observed  space  obiect  management  policy  which  could  ensure 
that  future  missions  would  be  protected  from  unacceptable  collision  hazards  in 
orbl  t. 
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APPENDIX 


SYNCHRONOUS  SATELLITE  CATALOG 
(28  April  1980) 
ALPHABETICAL 
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